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PREFACE 


Thil  projMt  mi  Initiated  bjr  the  U.  S.  Anor  Engineer  Oeodeiy, 
InteUlgenee  and  Mapping  Reaearoh  and  Derelopaant  kgmiey  (OIMRADA). 

It  is  being  pursued  under  Contract  0A'-ltli-009^EIf0-3767f  "Research  in 
Surveying,  Mapping  and  Qeodeay,"  authorised  by  and  negotiated  under 
Title  10,  U.  S.  C.,  Section  230U  (a)  (1)  and  Presidential  Proelanation 

29m. 


Starting  date  of  the  original  oontraot  waa  S^pteaber  U,  1958. 

The  work  reported  herein  waa  done  aainly  between  October  19ol  and 
November  19o2  under  subaequent  ■odifioations  (in  effect  Nuidbera  $ 
and  6)  of  that  oontraot.  Tendnal  date  of  Modification  No.  8, 
presently  in  force,  ia  June  30,  196b.  Execution  of  the  project  is 
under  the  direction  of  Mr.  Nathan  Flshel  of  OIMRADA,  Contracting 
Officer* a  Technical  Representative. 

This  project,  as  a  concern  of  the  Operations  Analysis  Laboratory 
of  the  Franklin  Institute  Laboratories,  has  been  under  the  supervision 
of  Mr.  Robert  S.  Orubaeyer,  Manager  of  that  Laboratory.  Mr.  Jaises  K.  Olela 
and  Dr.  John  E.  Merrill  carried  out  the  study  reported  here;  Dr.  WHliaa 
M.  Protheroe  served  as  consultant  on  various  aspects  of  the  problea. 

In  aooordance  with  the  pattern  established  in  Modification  No.  8 
of  the  basic  contract,  this  report  is  designated  Interia  Technical 
Report  1.  As  a  theoretical-ccaputational  investigation  of  the  relation 
between  the  aodel  chosen  to  represent  the  earth's  ataosphere  at  a 
given  tiae,  and  the  aaount  by  idiioh  an  inooalng  ray  of  light  would 
be  deviated  frea  its  original  direction,  it  is  closely  related  to 
the  study  of  the  19$9  ARDC  aiodel  presented  in  our  Quarterly  Report 
No.  6  and  that  of  aore  generaliaed  ataosphere  aodela  presented  in 
Part  I  of  our  Quarterly  No.  7;  for  this  reason  the  present  report  has 
III  aiqpended  to  the  newly-adored  sub-title. 
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SUMMARY 


This  report  continues  the  investigation  of  the  astronondeal 
refraction  in  spherically^etratified  model  atmospheres.  Previous^, 
the  study  of  the  refraction  in  the  1959  ARDC  model  was  presented 
Here  we  investigate  the  astronomical  refraction  in  tropical  and  arctic 
winter  atmosphere  models  at  zenith  distances  of  60^  and  85^*  The 
procedure  is  similar  to  that  used  in  the  study  of  the  1959  ARDC  model; 
the  results^  given  in  the  following  table,  show  satisfactory  agree¬ 
ment  with  their  respective  values  confuted  from  the  Willis  and  Pulkovo 
tables. 


Atmosphere 


Model 

z  =  60° 

z  =  85° 

Tropical 

93'.'oU3 

Bh2"292 

Arctic  Winter 

116. 

700.787 

Also  included  is  a  sample  calculation  for  a  homogeneous  atmosphere, 
considered  as  a  limiting  case.  Again  the  results  show  good  agreement 
with  Ptilkovo,  to  rather  considerable  zenith  distances. 

The  effect  of  injecting  water  vapor  into  the  dry  1959  ARDC  model 
was  examined.  The  procedure  consisted  in  revising  the  value  of  the 
index  of  refraction  at  each  level  where  necessary,  using  the  Barrell- 
Sears  expression  for  n.  The  results  indicate  that  for  a  zenith 
distance  of  60®  there  is  no  difference  in  the  contribution  to  the 
total  refraction  between  the  dry  and  humid  models  above  15  km  altitude, 
and  for  a  zenith  distance  of  65®,  there  is  no  difference  between  the 
two  models  above  30  km. 
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I.  INTRODUCTION 

This  report  discusses  work  on  the  prdblem  of  stnospherie  refraction 
beyond  that  reported  in  Quarterly  Report,  Q-A2217-6^^.  The  Min  por¬ 
tion  deals  with  the  investigation  of  the  amount  of  refraction  caused 
by  atmospheres  based  upon  tropical  and  arctic  winter  models;  a  simpli¬ 
fied  calculation  conpares  the  refraction  produced  by  a  homogeneous 
atmosphere  with  that  based  on  a  realistic  model;  and  a  preliminary 

discussion  of  the  change  in  refraction  caused  by  introducing  water  vapor 

[2] 

into  the  ARDC  model'-  -*  is  also  included. 

The  data  for  the  tropical  and  arctic  winter  (or  polar)  models  were 
obtained  from  tables  in  'Tlilitary  Standard  210  A,  Climatic  ExtrsMS 
for  Military  Equipment,"  published  by  the  DepartMnt  of  Defense, 

2  August  1957*  The  important  features  of  these  tables  are: 

1.  The  data  are  given  to  an  altitude  of  100,000  feet 

2.  Sea-level  tenqserature  is  +  32.1Ilt°C  for  the  tropical  model 
and  -  27°C  for  the  arctic  winter  atmosphere 

3*  Sea-level  pressure  is  760  nn  Hg  and  768.8062  mm  Hg  for  the 
two  models,  respectively 

U,  The  value  of  gravity  is  asau^  constant  throughout  both 
atmospheres,  (980. 61^  em/8»e^),  and 

5.  The  coqpoaition  of  both  atmospheres  is  assumed  constant 
throughout . 

The  tropical  atmosphere  represents  an  annual  Man  day  between  0°  and  20°  N 
latitude,  while  the  arctic  winter  profile  is  based  on  radiosonde  infonna- 
tion  obtained  at  seven  North  American  stations  north  of  latitude  60°  N. 


■1' 
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For  our  purposes  both  stnospheres  wsrs  sssunsd  to  be  vslld  st  the 

o’  " 

Istitude  for  vhioh  the  ARDC  stmosphere  use  constructed,  1(5  32  33  N. 

The  procedure  for  Integrstion  wss  simllsr  to  thst  used  for  the 
ARDC  model,  and  is  discussed  in  detail  in  Q-A2217-6  and  in  Section 
II.  A.  below.  The  main  problem  with  these  atmospheres  is  thst  the 
tabulations  extend,  for  perfectly  good  meteorological  reasons,  to  an 
altitude  of  only  30  km.  This  necessitated  the  grafting  of  an  ARDC 
top  onto  each,  in  order  to  Insure  that  all  the  atmospheric  layers  which 
might  contribute  to  the  refraction  were  Included  in  the  integration. 

The  technique  for  doing  this  is  also  given  in  Section  II.  A.  of  this 
report. 

Coogwrisons  with  values  from  the  Willis  ^  and  Pulkovo  tables 
were  made  and  the  results  and  analyses  of  these  are  included  in  a  later 
section. 


-2- 
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II.  INTEGRATION  OP  THE  TROPICAL  ATMOSPHERE  MODEL 

A.  Meehanlos  of  the  Integration 

The  procedure  for  the  integration  of  the  tropical  atmosphere  was 
similar  to  that  previously  used  in  the  ARDC  case.  The  basic  equation 
to  be  integrated  is  (from  Q-A2217-6,  page  3) 


The  value  of  "a"  was  taken  as  20890722.13  feet  and  the  zenith  angles 
for  which  the  ref ractlons  were  computed  were  60*^  and  6$°,  the  same  as 
for  the  ARDC  model.  N  (the  ratio  of  the  index  of  refraction  at  any 
height  to  the  index  of  refraction  at  the  surface)  was  calculated  for 
thirty- five  levels,  selected  at  equal  intervals  of  3000-foot  pressure 
heights  by  the  following  formula,  based  on  the  Barrell  and  Sears 
expression  for  n  at  any  temperature  and  pressure^ 

N^  =  0.9997382825  + 

where  t  =  teiqperature  (°C)  and  p  =  pressure  (mm  Hg)  obtained  from  the 
tropical  atmosphere  table  of  Military  Standard  210  A.  The  table 
actually  gives  p  in  inches  of  mercury,  but  this  was  converted  to  mm  Hg 
in  order  to  make  the  expression  for  N  coiqpatlble  with  that  previously  derived. 


O.38L6928729 


.p  +  (10U900  -  157t) 
(1  +  .00366lt)10^ 
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In  ordnr  to  instir*  that  no  sigtiifioant  amount  of  refraction  waa 
omitted  and  to  make  oonpariaona  with  other  nodela  more  convenient,  an 
extension  from  the  100,000-foot  top  of  the  tropical  atmosphere  table, 
to  90  km,  was  made.  This  was  done,  in  effect,  by  grafting  the  ARDC 
atmosphere  model  from  30.62?  km  to  90  km  onto  the  tropical  atmosphere. 

Two  assumptions  were  madet  (1)  as  before,  the  amount  contributed  to 
the  total  refraction  by  any  reasonable  atmospheric  model  above  90  km 
is  negligible  and  (2)  the  upper  part  of  the  ARDC  atmosphere  model  is  a 
very  close  approximation  to  what  the  upper  part  of  the  tropical  atmosphere 
model  would  be  if  it  existed.  This  procedure  requires  a  knowledge  of 
the  refraction  effect  at  the  Interface  where  the  two  models  were  Joined 
together. 

In  order  to  present  the  description  of  the  procedure  followed  as 
clearly  as  possible,  we  must  introduce  a  few  symbols  for  notation.  Let 
A  refer  to  the  ARDC  atmosphere,  which  was  reported  in  Quarterly  Report 
No.  6,  T  to  the  tropical  model,  TA  to  the  tropical  model  above  30  km, 

AW  to  the  arctic  winter  atmosphere  model,  and  AWA  to  the  arctic  winter 
model  above  30  km.  These  will  be  used  principally  as  subscripts; 
other  minor  notations  will  be  Introduced  as  necessary.  Also,  hence¬ 
forth  in  this  report  "Tropical  Atmosphere  Model"  will  be  taken  to 
mean  the  profile  based  on  the  data  from  the  Military  Standard  210  A 
tropical  atmosphere  table  up  to  30  km,  plus  a  top  based  on  the  1959 
ARDC  profile  from  30  km  to  90  km.  Similarly,  the  "Arctic  Winter 
Atmosphere  Model"  will  be  taken  to  mean  the  profile  based  on  the  data 
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from  the  Military  Standard  210  A  polar  atmoaphara  table  up  to  30  kn, 
plus  a  top  baaed  on  the  1959  ARDC  profile  from  30  kn  to  90  km. 

The  first  problem,  that  of  Integrating  the  ARDC  modal  above  30  km, 
using  the  tropical  sea-level  parameters  was  aooonplished  as  follows. 

In  Quarterly  Report  No.  6  the  integration  of  the  ARDC  atmosphere  was 
expressed  in  terms  of  the  variable  N^  ~  ^  index 

of  refraction  at  height  h,  and  n^^  is  the  index  of  refraction  at  the 
base  of  A.  occurs  not  only  in  the  numerator  (dN^),  but  also  Inside 
and  outside  the  brackets  in  the  denominator  in  equation  (1).  The 
integrand  of  this  expression  %ias  altered  for  the  present  case,  by 
putting 

*^OA 

N^^  =  X  ^  =  N^  X  1.0000156009  (3) 

and  (1  +  h/a)^  n|j^  =  (1  +  h/a)^  nJ  X  1.0000312022  (U) ' 

since  =  1.000277391 

and  =  1.000261786 

by  equations  (3,U)  of  Quarterly  Report  No.  6.  The  integration  for  the 
upper  part  of  the  tropical  atmosphere  was  actually  done  for  35  to 
90  km  only.  The  Interval  from  30.627  km  (100,000  ft)  to  35  km  is  a 
transition  Interval  treated  as  explained  below. 


-5- 
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The  8«oond  prolan,  refraction  by  the  tranaition  aone,  eaa  examined 
in  three  different  waya.  The  agreement  of  the  reaulta  indioatea  that 
any  one  of  the  three  ia  a  aatiafaotory  approach  to  our  problem. 

The  first,  and  moat  obvious  approach,  is  by  means  of  Snell's  law. 

In  Figure  1  let  i^  be  the  angle  that  the  outgoing  ray  makes  with  the 
vertical  at  the  transition  sons,  treated  as  a  thin  shell  or  interfaoe 
of  constant  index  n^.  Then  nQipr^  sin  a  =  n^r^  sin  1^.  Similarly,  if 
i^  is  the  angle  the  outgoing  ray  makes  with  the  vertical  at  the  inter¬ 
face,  assuming  an  index  n^  on  the  upper  aide  of  the  surface,  then 
n^^TQ  sin  z  s  n^r^  sin  i^  and  therefore  np  sin  i^  =  n^  sin  1^  which  is 
the  ordinary  law  of  refraction  at  an  Interface  between  media  with  indices 
of  refraction  n^  and  n^.  Solution  of  these  equations  results  in  a 
value  of  0.0131  for  i^  -  i^^  at  z  =  60®  and  0^0568  at  z  =  85°.  These 
quantities  represent  the  amount  by  idiich  the  outgoing  ray  refracts 
towards  the  normal  in  crossing  the  Interface. 


-6- 
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The  second  nethod  egein  considers  the  trensition  sone  ss  en 
extreswly  sbibII  interval,  that  is,  an  interface,  the  lower  surface  of 
which  is  the  top  of  the  tropical,  and  the  upper  surface,  the  bottom 
of  the  TA  atmospheric  model.  This  in  effect  puts  -  h^  and  hence 
p?  =  Pj^  since  the  tropical  atmosphere  table  is  given  with  ICAO  pressure 
height  as  an  argument.  However,  and  therefore  n^  /  n^  and 

N.  N„.  In  fact,  N-.  =  -iii  =  0.9997U19033  and  N.„  =  ~  =  0.9997Ul9i*12. 

A  i  1  Hqy  Hi 

tl 

This  procedure  results  in  refractions  in  the  amount  of  0*0133  for 

o  tt  . 

2  =  60  and  0.0602  for  z  =  85  .  Note  further  that  since  ^n^^ 
the  ray  traced  outwards  bends  toward  the  vertical  in  crossing  the 
interface^  as  previously  concluded.  Thus,  in  effect,  a  negative 
increment  is  Introduced  into  the  sumnatlon  of  the  refraction  integral* 
The  third  procedure  Joins  the  atmospheres  at  a  point  where  the 
tonperatures  come  together*  This  involves  an  extrapolation  of  the 
temperature  gradient  of  the  tropical  atmosphere  and  asstimes  that  both 
gradients  are  linear  when  expressed  as  a  function  of  geopotential 
feet.  This  is  not  quite  true,  but  the  deviation  is  negligible  for  our 
purposes*  The  presstires  at  this  altitude  will  also  be  equal,  hence 
the  n*s  and  N^s*  With  these  asaui^tions  the  temperatures  were  found 
to  come  together  at  3U»858  geopotential  km  (or  35*05  geometric  km)* 

The  coonon  temperature  at  this  level  is  --  26*93^0  and  the  ICAO  pressure 

is  U*3U8  mm  Hg*  An  N  for  this  level  was  calculated  and  the  computation 

a 

of  the  refraction  for  the  30.62?  to  35*050  kra,  35*050  to  1*0  km,  and  the 
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UO  to  90  kxn  intervals  constitutes  the  third  method  of  grafting  the 
two  atmospheres  together*  The  results  for  this  procedure^  for  both 
the  60^  and  the  85^  zenith  angle  cases^  are  given  in  Section  B* 

B*  Results 

The  refraction  produced  by  an  atmosphere  of  the  model  assumed^ 
and  confuted  in  the  manner  described  in  Section  II.A*  is  given  in 
Table  1  for  the  two  zenith  distances,  60®  and  85^*  The  left  side  of 
the  table  presents  the  refraction,  accumulated  from  the  effective  top 
of  the  TA  model,  down  to  any  level,  h,  terminating  at  the  upper  side  of 
the  interface*  The  right  side  of  the  table  presents  the  refraction, 
accumulated  from  the  top  of  the  tropical  atmosphere  proper  down  to  any 
level  h#  The  entries  in  the  interval  column  were  obtained  by  converting 
the  entries,  taken  at  3000-foot  intervals  in  the  altitude  coliamn  of  the 
tropical  atmosphere  table  into  kilometers  (1  ft.  *  .30118  meters  exactly, 
ICAO  Extension,  1959)#  Also  included  in  this  table  are  values  of  the 
refraction  at  the  interface  calculated  by  the  third  method  described 
previously.  Table  2  summarizes  the  results  for  the  integration  of  the 
tropical  atmosphere. 
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Table  1 


REFRACTION  ACCUMULATED  DOWN  TO  ANY  LEVEL,  h, 
IN  THE  TROPICAL  ATMOSPHERE  MDDBL 


Level 

Level 

h  (km) 

z  ■  60® 

z  -  85° 

h  (km) 

z  ■  60® 

z  •  85^ 

90-85 

H 

11 

It 

n 

0.0003 

0.0009 

30.7-30.1*  0.0601* 

0.2736 

80 

.0013 

.001*3 

29.5 

.2711 

1.2336 

75 

.003U 

.0111* 

28.5 

.5110 

2.3362 

70 

.0075 

.0257 

27.6 

.7969 

3.6607 

65 

.011i7 

.0513 

26.7 

1.1321* 

5.2293 

60 

.0271 

.0966 

25.7 

1.5326 

7.U65 

55 

.01*71 

.1717 

21*.8 

1.9981 

9.3320 

50 

.081*2 

.311*8 

23.9 

2.51*71* 

11.9703 

li5 

.1582 

.6099 

23.0 

3.1802 

15.0381 

Uo 

.31U3 

1.2539 

22.1 

3.9266 

18.6906 

35 

.6509 

2.6953 

21.2 

U.8107 

23.0580 

30.627 

1.2783 

5.1*6U9 

20.3 

5.8895 

28.1*380 

19.1* 

7.1571 

31*.  8202 

Uo-35.05 

0.3319 

1.1*208 

18.5 

8.61*7U 

1*2.391*7 

J.05-30.627 

0.6188 

2.7U99 

17.7 

10.1*0Ul 

5l.li076 

16.9 

12.1*691* 

62.101*2 

16.1 

11*.261*2 

71.1*882 

15.2 

16.1856 

81.6337 

ll*.3 

18.31*61 

93.1598 

13.1i 

20.7690 

106.2266 

12.5 

23.1*809 

121.0185 

11.6 

26.51*52 

137.9331* 

10.6 

29.9191 

156.7897 

9.6 

33.6107 

177.6871 

8.7 

37.6293 

200.7345 

7.7 

1*1.9815 

226.0315 

6.8 

1*6.7091* 

253.8901* 

5.8 

51.8293 

281*.U81*6 

1*.8 

57.3296 

317.8255 

3.9 

63.2851 

351t.U587 

2.9 

69.6666 

394.3059 

1.9 

76.5239 

437.7861 

1.0 

83.9065 

485.3395 

0.0 

91.7782 

536.8664 
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SUItfARY  OF  RESULTS,  NUMERICAL  INTEGRATION 
OF  THE  TROPICAL  ATMOSPHERE  MOIEL 


z  •  60° 


Level  (kmj 

Hal 

a 

bhod 

b 

Level  (km) 

Method 

0 

90-30.627 

Interface 

30.627-0 

l!l2783 
-  0.0131 
91.7782 

l!l2783 
-  0.0133 
91.7782 

90  -  Uo 
^  lO  -  35.05 

i  mh : 

o"311t3 

mi 

_ 91*22a2_ 

Totals  93"Olt3U  93"oli32 

Total  93"o1*32 

■  -  650 


Level  (tan) 

Method 

a  b 

Level  (km) 

Method 

e 

90-30.627 

Interface 

30.627-0 

5“U8U9 
-  0.0588 
536e866U 

5'.'48U9 
-  0.0602 
536.866U 

90  -  liO 

UO  -  35.05 
35.05  -  30.627 
,  30.627  -  0 

l"2539 
2.7U99 
1.U208 
_ S36.886U 

Totals  5U21I2925  5U2ll29U 

Total  5U2II29IO 

Method  a  -  Snell's  law 

Ifethod  b  -  Interface  as  a  thin  layer 

Method  c  -  Extrapolation  of  temperature  gradient 


C*  Comparisons 

Table  3  shows  the  con^jarisons  of  the  refraction  obtained  by  the 
integration  of  this  model  atmosphere,  with  that  obtained  fron  the  Willis 
and  PxilkoTo  systems,  using  the  same  values  for  the  parameters  at  the  « 
base.  A  discussion  of  these  results  will  be  postponed  until  the  arctic 
winter  data  are  added  in  a  later  section  of  the  report. 

Table  3 

COMPARISONS  OF  THE  REFRACTIONSt  ‘nUDPICAL  INTEGRATION^  WILLIS, AND  PULKOVO 


Z  ■  60® 

z  -  85° 

Tropical 

93l!oU3 

5U27292 

WiUis 

92.990 

5U2.678 

Pulkovo 

92.9U6 

5U2.062 
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III.  ZNTECSUTION  OF  THE  ARCTIC  WINTER  ATMOSPHERE  MOIEL 
A.  Mschsnies  of  the  Integration 

The  procedure  for  the  integration  of  the  arctic  winter  atmospheric 
model  was  the  same  as  that  for  the  tropical  model*  The  8ea«-level  values 
of  the  temperatures  and  pressure  are  -27®C  and  768.8672  mm  Hgj  respectively* 
The  AW  atmosphere  was  assumed  to  be  valid  for  the  latitude  of  the  ARDC 
roodelj  hence  the  value  of  the  earth’s  radius  was  unchanged*  The  value 
of  la000328653>  was  confuted  from  equation  (h)  of  Quarterly 

Report  No*  6*  The  were  confuted  from  the  following  expression! 

-  o.9996nusso  -  2-3awn?ei  Cp  .  juph?,  -  (5, 

(1  +  •00366lt)10^ 

An  ARDC  top  was  oon9)uted  for  the  AW  models  in  much  the  same  way  as 
for  the  tropical  model*  In  order  to  equalize  the  pressures  in  the  two 
atmospheres^  it  was  necessary  to  lower  the  ARDC  levels^  in  effect^  by 

I 

1*3  Rm  resulting  in  a  new  set  of  N  a, 

The  Interface  of  this  ecaibined  atmospheric  model  is  at  29.337  geometric 
kilometers,  or  96,2U9  geometric  feet  (as  given  by  the  extreme  right-hand 
column  of  the  ajrctic  winter  atmosphere  table).  The  refraction  at  this 
level  was  computed  by  the  first  two  methods  previously  used,  namely  by 
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Snell's  law  and  by  the  introduction  of  a  smallf  thin  layer.  The  first 
procedure  resulted  in  reAractions  of  o!!lUl2  and  o!!6L32  for  z  •  and 
z  ■  85^  respeotivelyt  while  the  second  method  resulted  in  refractions 
of  o!llill  and  o!l61i90  for  z  -  60°  and  z  -  85°  respectively,  in  all  cases^ 
the  ray  traced  outwards  from  the  observer  bends  away  from  the  normal 
in  crossing  the  interface,  opposite  to  the  situation  found  for  the 
tropical-ARDC  model. 

B.  Results 

The  refraction  produced  by  an  atmosphere  of  the  model  assumed,  and 
conf>uted  in  the  manner  described  in  Section  III.A.,  is  given  in  Table  U 
for  the  two  zenith  distances,  60°  and  85°.  The  left  side  of  the  table 
presents  the  refraction,  accumulated  from  the  effective  top  of  the  AWA 
model,  down  to  any  level,  h,  terminating  at  the  upper  side  of  the  inter¬ 
face.  The  right  side  of  the  table  presents  the  refraction,  accumulated 
from  the  top  of  the  arctic  winter  atmosphere  proper  down  to  any  level, 
h.  Also  included  in  this  table  is  the  value  of  the  refraction  at  the 
interface  calculated  according  to  the  second  method  described  above. 
Table  5  summarizes  the  results  for  the  integration  of  the  arctic  winter 
atmosphere. 
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Table  U 

REFRACTION,  ACCUMUUTED  DOWN  TO  ANY  LEVEL,  h, 
IN  THE  ARCTIC  WINTER  ATMOSPHERE  MODEL 


Level 

Level 

h  (km) 

z  -  60° 

z  -  85° 

h  (km) 

z  ■  60® 

z  ■  85° 

88.7-83.7 

o!o003 

o"ooii 

29.3-29.1  o!!o650 

o!!2994 

78.7 

.oolU 

.0045 

28.2 

.2797 

1.3004 

73.7 

.0034 

.0115 

27.4 

.5294 

2.4652 

68.7 

.0075 

.0259 

26.5 

.8215 

3.8386 

63.7 

.0147 

.0518 

25.7 

1.1612 

5.4494 

58.7 

.0272 

.0975 

24.8 

1.5532 

7.3238 

53.7 

.0472 

.1734 

23.9 

2.0095 

9.5249 

1»8.7 

.0843 

.3183 

23.0 

2.5362 

12.0888 

U3.7 

.1584 

.6171 

22.1 

3.1414 

15.0619 

38.7 

.3147 

1.2700 

21.2 

3.8427 

18.5396 

33.7 

.6516 

2.7324 

20.3 

4.6500 

22-.  5820 

29.337 

1.2797 

5.5669 

19.4 

5.5838 

27.3039 

18.5 

6.6554 

32.7774 

Interface 

♦0.1411 

♦0.6490 

17.6 

7.8956 

39.1769 

16.7 

9.3262 

46.6374 

15.8 

10.9720 

55.3120 

14.9 

12.8774 

65.4651 

14.0 

15.0723 

77.2926 

13.1 

17.5986 

91.0615 

12.2 

20.5195 

107.1669 

11.3 

23.8910 

125.9776 

10.4 

27.7664 

147.8621 

9.5 

32.1823 

173.1035 

8.6 

37.0858 

201.4762 

7.7 

41.7170 

228.6024 

6.8 

46.7648 

258.5358 

6.0 

52.2772 

291.6337 

5.1 

58.2356 

327.8642 

4.3 

64.7189 

367.7926 

3.4 

71.7291 

411.5262 

2.5 

79.6504 

461.5899 

1.7 

89.1921 

522.6707 

0.9 

99.8071 

591.4708 

0.1 

U3.9409 

684.1398 

0.0 

n5.5229 

694.5439 
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Table  5 


SIMIARY  OF  RESULT^  NUMERICAL  INTEORATION 
OF  THE  ARCTIC  WINTER  ATMOSPHERE  MODEX 


z  -  60° 


Level  (km) 

Method 

Method 

a 

b  . 

87.7  -  29.337 

l'.'2797 

Interface 

0.11*12 

Total  Il6.9l;38  116.91*37 


z  -  85® 


Level  (km) 

Method 

Method 

a 

b 

87.7  -  29.337 

5'.'5669 

5'.’5669 

Interface 

+  0.61*32 

+  0.61*90 

29.337  -  0 

69U.5739 

69l*.5739 

Total 

700'.' 781*0 

70o!7878 

Method  a  -  Snell's  law 

Ifethod  b  -  Interface  as  a  thin  layer 


C.  Con^jarisons  and  Analyses 

Table  6  and  Figures  2  and  3  svunmarize  the  refraction  confutations. 
In  Table  6  the  integrations  of  the  three  model  atmospheres  are  confared 

a 

with  the  corresponding  values  of  the  refractions  obtained  by  the  Willis 
and  Pulkovo  methods*  To  repeat ^  the  same  basic  parameters  were  used  in 
the  Willis  and  Pulkovo  methods  as  in  the  numerical  integrations  of  the 
atmospheric  models^  namelyt 

1*  Wavelength  of  lights  *5753  microns 

2*  Standard  gravity  at  latitude  ii5^32*33'*>  980.665  cm  sec"^ 

3*  No  water  vapor  in  any  of  the  models*  In  addition^  the  surface 
conditions  for  each  atmosphere  were  as  follows  t 
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■15' 


h  (km) 
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Atmospheric  Model 


T(°C)  p  (mm  Hg) 


I-A2217-i 

(Hq  -  1)  X  10^ 


Tropical  32*llU 

ARDC  +  15 

Arctic  Winter  -  27 


760 

760 

768,8072 


261.786 

277.391 

328.653 


At  z  -  60°,  the  results  are  fairly  consistent.  The  integration  method 

n 

in  all  three  cases  gives  about  0.05  more  refraction  than  Willis*  The 

tt 

integration  also  gives  about  0*1  more  refraction  than  Pulkovo^  except  in 

ft 

the  ajxtic  winter  case  where  it  gives  about  0*16  more*  Willis  gives 
about  o!!o5  more  refraction  than  Pulkovo  for  the  first  two  model  atmospheres. 

It 

and  ‘t*  0*11  more  for  the  arctic  winter  atmosphere*  One  might  place  the 
blame  for  the  deviation  in  the  cold  atmosphere  case  on  the  Pulkovo  tables, 
but  this  is  difficult  to  Justify  in  view  of  the  fact  that  these  tables 
are  observationally  based|  their  standard  conditions  being  representative 
of  a  high  northern  latitude*  At  a  zenith  angle  of  85^  the  numerical 
values  are  less  well  stabilized  because  any  error  is  multiplied  by  a 
factor  of  almost  ten  (with  respect  to  the  data  for  z  «  60^)  and  because 
we  are  working  at  extreme  limits  in  both  tables,  where  even  rounded-off 
decimals  may  be  significant*  For  exair^^le,  within  the  Willis  tables,  an 
uncertainty  of  ~  0?008  can  be  e3q>ected  in  the  value  of  the  refraction 
at  a  zenith  distance  of  85^9  <^e  primarily  to  the  restriction  of  the 
number  of  decimals  given  for  computing  two  of  the  basic  parameters,  f 
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and  log  n  •  An  uncertainty  of  5  x  10~^  In  the  value  of  log  R  derived 
0  0 

from  the  Pulkovo  tablee^  gives  an  tmeertalnty  of  approximately  the 
same  order  of  magnitude  in  the  refraction  as  that  found  for  the  Willis 
tables* 

Table  6 


SUMMARY  OF  REFRACTION  COMPUTATIONSi 
THE  INTEGRATIONS  OF  THE  THREE  MODEL  ATMOSPHERES 
AND  THE  CORRESPONDING  WILLIS  AND  PULKOVA  VALUES 


Method 


Tropical 


ARDC 


60“ 


Arctic  Winter 


Integration 

93.0U3 

98!!620 

116.91*1* 

Willis 

92.990 

98.563 

116.890 

Pulkovo 

92.9U6 

98.512 

116.781* 

tl 

II 

II 

I  -  W 

♦  .053 

+  .057 

+  .051* 

I  -  P 

+  .097 

+  .108 

+  .160 

W  -  P 

+  .Ohh 

+  .051 

+  .106 

z 

O 

CO 

■ 

Integration 

Sh2'.202 

578.811* 

700! 787 

Willis 

5U2.682 

579.026 

698.888 

Pulkovo 

5U2.138 

578.512 

698.933 

n 

N 

n 

I  -  W 

-  .390 

-  .212 

+  1.900 

I  -  P 

+  .13U 

+  .302 

+  1.855 

W  -  P 

+  .5Ui 

+  .511* 

-  .01*5 

At  a  zenith  distance  of  85°»  Table  6  shows  that  the  refraction 
obtained  by  integrating  both  the  tropical  and  ARDC  model  atmospheres 
was  less  than  that  obtained  by  the  Willis  method  and  more  than  that 


obtained  from  the  Pulkovo  tables.  The  results^  though  not  in  as  good 
agreement  as  those  for  a  zenith  distance  of  6o“}  are  consistent  when 
allowance  is  made  for  the  fact  that  at  such  a  large  zenith  distance^ 
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the  effects  of  the  small  differences  in  the  atmospheric  models  (those 
on  which  the  tables  are  based)  become  highly  magnified*  Willis  gives 

w 

about  0*5  more  refraction  than  Pulkovo  when  the  surface  conditiona 
are  those  of  the  tropical  and  ARDC  atmosphere  models*  However^  for 
the  arctic  winter  atmosphere  conditions^  the  integration  gives  nearly 
two  seconds  more  refraction  than  either  the  Pulkovo  or  Willis  methods^ 
and  the  latter  two  come  more  into  agreement,  Pulkovo  giving  about  0*05 
more  refraction  than  Willis* 

In  order  to  see  how  our  three  cases,  two  of  which  lie  near  the 
extremes  in  their  sea  level  temperature  parameters,  fit  into  the  over¬ 
all  pattern  of  atmospheric  refraction,  a  study  was  made  of  the  refraction 
for  the  range  of  base  conditions  likely  to  be  encountered  on  the  earth's 
surface*  The  basic  qixantities  on  tdiich  these  refraction  confutations 
were  made  aret 


t  -  +30,  +15,  0,  -15,  -30  (°c) 

p  -  720,  760,  780  (mm  Hg) 

z  -  60°,  70°,  80°,  85° 

Humldityt  0% 

The  results,  shown  as  differences  between  the  Willis  and  Pulkovo  values 
of  the  refraction  (in  the  sense  VHllis  minus  Pulkovo)  are  given  in 
Table  7  and  are  plotted  in  Figure  U  as  functions  of  the  teiq)erature  for 
each  zenith  angle*  The  run  of  the  differences  forms  a  consistent 
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progression  from  z  ■  60°  to  z  ■  8$°*  However,  the  large  range  of  the 

differences  encountered  at  the  zenith  distance  of  85^,  makes  It  fairly 

clear  that  at  large  zenith  distances  one  should  realize  that  In  speak- 

N  n 

Ing  of  orders  of  accuracy  such  as  0.001  or  even  0.01,  one  must  be  very 
careful  In  making  calculations  with,  and  Interpretations  of,  the  numbers 
which  are  obtained  from  any  existing  method  of  refraction  oonyiutatlon* 
New  model  atmospheres  have  been  recently,  or  are  about  to  be, 
published.  Any  further  study  of  refraction  along  the  lines  pursued 
here  should  naturally  make  use  of  them.  Because  of  time  limitations 
brought  about  by  the  pressure  of  other  problems,  there  still  remain 
some  closely-related  questions  to  be  answered  -  for  exan^le,  the  effect 
of  refraction  at  still  larger  zenith  distances  where  the  basic  Integral 
converges  even  more  slowly  for  any  realistic  model. 

Table  7 


DIFFEEIENCES  IN  REmCTION  BASED  ON  THE  WILLIS 
AND  PULKOVO  SYSTEMS,  IN  THE  SENSE,  W-P 


t 

oc 

p 

Bim  Hg 

z  «  60® 

z  -  70° 

O 

O 

CO 

■ 

z  -  85* 

-30 

720 

♦ro89 

♦.15U 

+.161 

-!i80 

760 

♦.093 

+.158 

+.178 

-.143 

780 

♦.095 

+.16U 

+.183 

-.128 

-15 

720 

♦.073 

♦.13U 

+.193 

♦.178 

760 

♦.077 

+.131 

+.199 

♦.188 

780 

♦.079 

+.138 

+.209 

+.197 

0 

720 

♦.062 

+.110 

+.196 

♦.U30 

760 

♦.065 

+.1U 

+.200 

+.10.8 

780 

♦.065 

+.115 

+.210 

+.419 

♦15 

720 

♦.0U8 

+.087 

+.179 

+.540 

760 

♦.051 

+.085 

+.181 

♦.514 

780 

♦  .052 

+.090 

+.18U 

♦.503 

♦30 

720 

♦.037 

+.067 

♦.139 

+.600 

760 

♦.037 

+.062 

♦.133 

♦.549 

780 

♦.038 

♦,06U 

♦.13U 

♦.528 

-21- 


THE  FRANKLIN  INSTITUTE  •  labmtoiiniwJiimanAian^ 

I-A  2217-1 


IV.  REFRACTIOM  BY  A  HOMOQENEOtS  ATHDSPHERE 


A  very  eiiqile  oalouletlon  wee  made  to  deteraina  the  magnitude  of 
the  refraction  caueed  by  a  haaogeneoua  atmosi^ere  thought  of  as  a 
limiting  ease.  "Homogeneous*'  was  taken  to  mean  that  the  dependenoe  of 
tesperature*  pressure  and  density  on  height  was  such  that  the  refractive 
Index  remained  constant  throughout.  The  height  of  such  an  atmosphere ^ 
with  a  sea-level  teoperature  of  Is  given  by  Allen  [^]  as  8.U30  km. 

A  series  of  calculations  were  made  for  varloius  zenith  angles  based  on 
the  following  formula  derived  by  use  of  Snell's  lawt 


sin"^ 


Ug  a  sin  z 
a  4-  h 


-  sin 


-1 


a  sin  z 
a  +  h 


(7) 


tdiere,  as  before,  n  ■  Index  of  refraction  at  the  earth's  surface, 

1.000277391 

a  -  radius  of  earth,  6367.U9190  km 

z  ■  angle  between  the  zenith  and  the  direction  of 
the  ray  at, the  observer 

h  ■  altitude,  or  In  this  case,  the  height  of  the 
homogeneotu  atmosphere,  8.U30  km 

The  results  are  shown  in  Table  8,  together  with  corresponding  values 
coi^puted  from  the  Rilkoeo  tables.  It  appears,  rather  surprisingly,  that 
down  to  an  apparent  zenith  distance  of  80**,  even  the  extremely  artificial 
hoaMgeneous  atmosidiere  provides  a  value  of  the  reftractlon  actually  good 
enougi  for  many  ordinary  purposes  such  as  reduction  of  survey  plates, 
and  so  forth. 
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Table  8 

REFRACTION  BT  A  HOMOCffiNEOUS  ATMOSPHERE 


Apparent 

z.  d. 

Honog 

V*lk 

H  -  P 

o 

M 

n 

M 

0° 

0 

0 

0 

10 

10.076 

10.063 

+  0,013 

20 

20.79U 

20.770 

O.O2I* 

30 

32.977 

32.91*2 

0.035 

1*0 

U7.906 

U7.858 

O.OU8 

50 

67.985 

• 

67.910 

0.075 

60 

98.621 

98.512 

0.109 

70 

155.610 

155,1*1*6 

0.161* 

75 

209,821* 

209.61(8 

+  0.176 

80 

312.397 

312^1*79 

-  0.082 

85 

570.897 

578.512 

-  7.615 

23- 


THE  FRANKLIN  INSTITUTE  •  UAmttnin  for  RtnarA  and 

I-A2a7-1 


V.  THE  EFFECT  OF  MOISTURE 

The  AROC,  tropical,  and  arctic  winter  model  atmospheres  are  based 
upon  the  assunption  of  dry  air.  It  is  well  known  that  humidity  affects 
the  astronomical  refraction;  both  the  Pulkovo  and  Willis  methods  of 
computation  make  allowance  for  the  influence  of  moisture.  This  allowance 
is,  however,  based  upon  the  moisture  content  of  the  air  in  the  immediate 
vicinity  of  the  observer.  To  test  the  correction  for  water  vapor  it 
was  decided  to  introduce  moisture  into  the  ARIX!  model  atmosphere  at  all 
levels  and  to  re-integrate  the  refraction;  this  was  then  to  be  compared 
with  the  values  predicted  by  the  Willis  and  Pulkovo  calcxilations.  It 
was  assumed  that  the  moisture  profile  given  by  Murcray  et  al j  Is  re¬ 
presentative  of  the  normal  atmospheric  conditions  when  moisture  is 
present.  Since  large  variations  in  the  moisttire  content  of  the  air  with¬ 
in  a  few  kilometers  of  the  surface  can  exist,  it  was  decided,  for  ease 
of  computation,  to  use  a  profile  consisting  of  four  straight  line  seg¬ 
ments  to  represent  the  mixing  ratio  of  water  to  air  by  mass  as  a  function 
of  height.  The  assumed  function  and  the  original  Murcray  profile  are 
shown  in  Figure  5, 

The  introduction  of  water  vapor  into  the  ARDC  model  entails  a  small 
adjustment  in  one  or  more  of  the  three  parameters^  ten^erature,  pressure, 
and  density,  describing  the  a^aosphere  as  a  function  of  heic^t.  This 
is  caused  by  the  modificatim  in  the  mean  molecular  weight  of  the  air 
brought  about  by  the  admixture  of  water.  Since  any  change  in  density 
would  in  turn  require  a  change  in  the  pressure  profile  of  the  model. 
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It  was  dacldsd  to  introduoa  the  amount  of  water  vapor  required  by  the 
asstoMd  nixing  ratio  profile  in  such  a  manner  that  the  density  of  the 
moist  atmosphere  would  resialn  the  sane  as  that  tabulated  in  the  dry  ARDC 
model.  This  procedure,  which  retains  the  pressure  profile  of  the  dry 
AR1)C  and  concentrates  the  adjustment  into  the  temperature  profile, 
seemed  to  be  the  one  which  would  be  likeliest  to  bring  the  tabulated 
ARDC  model  close  to  the  reality  of  a  humid  atmosphere. 

If  the  air  is  taken  to  obey  the  perfect  gas  law  in  accordance  with 
the  assumption  used  to  construct  the  ARDC  model,  then  the  absolute  tenqpera- 
ture  of  the  air  at  any  level  is  given  by 

T-J|m  (8) 

where  T  ■  tenperatiure  in  °K 

P  «  the  pressure  in  mm  Hg 
p  «  the  density  in  grams/cc 

K  -  the  perfect  gas  constant  in  appropriate  units 
N  «  the  mean  molecular  weight  (dimensionless) 

Ft*on  the  definition  of  mean  molecular  weight  and  the  mixing  ratio  It 
is  possible  to  write  for  each  level 

M  -  28.966  -  17.520W  (9) 

where  w  *  the  mixing  ratio  in  grams  of  H2O  per  gram  of  air  and  28.966  is 
the  mean  molecular  wei^t  of  dry  air  for  the  coaposition  assumed  in  the 
ARDC  model. 
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The  perfect  gee  lew  then  gives 


T  -  J  5  (28.966  -  17.520W) 


(W) 


The  tenqperature  of  the  dry  AROC  atmosphere  at  any  level  can  be  written 
as 


T 


d 


ifd 

|[5^"d 


(11) 


where  the  subscript  d  indicates  the  value  of  the  quantity  designated  at 
the  given  level  in  the  1959  (dry)  ARDC  tables*  Taking  the  ratio  of  the 
two  expressions  and  noting  that  P  ■  and  P  ■  by  the  assuzi^tion 
above  and  that  «  26.966^  the  temperature  of  the  moist  air  may  be 
written  in  terms  of  the  "dry  ARDC"  tenperature  and  the  mixing  ratio  as 


T  •  (1  .  0.60l485w)  (12) 

• 

This  gives  the  value  which  the  temperature  of  the  moist  air  must  assume 
at  each  level  in  order  that  its  density  and  therefore  its  pressure  at 
that  level  shall  be  the  sane  as  in  the  ARDC  tables. 

With  the  neWf  or  modified  teii|>erature«  it  is  necessary  to  recompute 
the  refractive  index  at  each  level  in  the  atmospherej  the  new  value  tms 
taken  to  be  that  given  by  the  ejqpression  of  Barrell  and  Sears  for 
the  refractive  index  of  moist  air,  namely 


(n  -  1)  -  ^nP  *  ^tP)  _  f  (1  ♦  2^^)  %  -  K« 
1  ♦  at  1  ♦  ot 


(13) 
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Th«  quantities  in  this  eiqpresslon  are  defined  asi 

-  [0.378125  +  0.0021laltX  ♦  0.00001793X  “^]  x  10“^ 

K'  -  [0.3159  ♦  0.002963X  ]  X  10“^ 

p  ■  total  pressure  in  nm  Hg 
a  ■  constant  ■  3*661  x  10“^ 
jS^  -  (1.0it9  -  0.0l57t)  X  10"^ 

t  -  ten^erature^  ^Celsius 

The  wavelength  of  the  lic^t  was  taken  to  be  X  -  0.5753  microns  as  in  our 
other  studies  of  atmospheric  refraction,  giving  the  values 

-  0.38U75875  x  10"^ 

and 

K*  -  0.32U852U7  X  10“^ 
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The  second  term  of  the  Barrel  and  Sears  eiqpression  represents  the 
direct  influence  of  the  water  vapor  upon  the  value  of  (n  -  !)•  The 
value  of  (n  >  1)  is  also  atodified  indirectly  at  each  level  of  the  humldi> 
fled  model  atmosphere,  due  to  the  modification  in  the  teiqperature 
described  above* 

While  it  is  possible  to  confute  (n  -  1)  directly  at  each  level 
from  the  Barrel-Sears  formula,  it  may  also  be  e:q>ressed  as  an  esqpansion 
of  this  formula  about  the  value  at  t  ■  t^,  w  >  f  ■  0*  This  permits 
additional  insight  into  the  effects  upon  the  index  due  to  the  two  causes, 
change  in  tenperature  and  inherent  influence  of  water  vapor  upon  re¬ 
fraction.  l^on  expansion  and  some  simplification  the  expression  for 
-  1)^  the  index  of  refraction  for  moist  air^  may  be  written  as 


0.62197  -  w 

iy:E2i97  ” 


K 

to 

f(%  ♦  k')( 

li.,  ‘  1 

"  1  .  atd  ^ 

J 

1 

(lU) 


where  ■  the  index  of  refraction  for  the  given  level  in  the  ARDC  model 
y  -  9.1i96l  wT^  X  10"^ 

A  -  2.211*356  wTjj  X  10“^ 

and  T^  ••  absolute  teii9)erature  of  the  ARDC  model  at  the  corresponding 
height.  This  expression  reduces  to  that  used  for  dry  air  in  our  previous 
study  M  and  is  coxrect  to  much  better  than  one  percent  for  all  reason¬ 
able  values  of  w,  t^  and  p  encountered  in  the  problssi* 
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An  analysis  of  this  expression  indicates  that  the  refractive 
index  for  the  moist  atmosphere  will  always  be  greater  then  that  of 
the  AROC  atmosphere  for  corresponding  heights.  The  values  of 
,  (n^  -  1)  were  computed  at  selected  heights  and  compared  with  the 
values  of  (n^  -  1)  already  available  at  those  heights.  The  results 
of  the  calculations  are  given  in  Table  9. 


Table  9 


COMPARISON  OF  INDICES  OF  REFRACTION  AT  SELECTED  LEVELS  IN  THE  MODIFIED 
(WET)  ARDC  ATMOSPHERE  AND  IN  THE  1959  (l»Y)  ARDC  ATMOSPHERE 


Wtl 

AHDC 

Dry  ARDC 

Difference 

h 

w  X  10^ 

(n  -  1 )  X  10^ 

'  w 

(n^  -  1)  X  10^ 

(\  - 

55 

0.100 

138,01 

138 

0.01 

50 

0.100 

2li5.01 

21*5 

0.01 

30 

0.100 

U,0U2.11* 

l*,0l*2 

O.ll* 

25 

0.032 

9,198.10 

9,198 

0.10 

20 

0.010 

20,122.07 

20,122 

0.07 

15 

0.010 

Uh, 080.07 

1*1*,  080 

0.07 

10 

0.028 

93,610.9 

93.610 

0.9 

9 

0.0U6 

105,739.7 

105,738 

1.7 

8 

0.077 

119,039.2 

119,036 

3.2 

7 

0.129 

133,589.1 

133,583 

6.1 

6 

0.215 

lh9,U69.U 

11*9.U58 

U.l* 

5 

0.359 

166,760.2 

166,739 

21.2 

U 

0.599 

185,560.2 

185,521 

39.2 

3 

1.000 

205,957.9 

205,885 

72.9 

2 

1.670 

228,060.9 

227,926 

13U.9 

1 

2.780 

251,975.3 

251,727 

21*8,3 

0 

Ua6U0 

277,81*8.8 

277,391 

1*57.8 

It  is  obvious  from  an  Inspection  of  the  table  that  most  of  the 
change  in  refraction  resulting  from  introduction  of  the  water  vapor 


occurs  close  to  the  surface  (more  than  70%  of  the  change  is  within  the 
lower  2  km)  and  that  practically  all  significant  change  has  occured 
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witMn  the  firat  half  dosan  or  so  kilonetera*  This  nay  be  taken  aa  a 
reaaonable  esqplanatlon  of  the  known  fact  that  the  "humidity  oorreotlona" 
uaed  with  the  Pulkovo  and  WllUa  tablea  (though  baaed  only  on  the  local 
aurfaee  humidity)  are  valid  for  moat  purpoaea*  It  is  of  course  possible 
that  a  suitably,  but  slightly  augmented  value  of  n  based  i;^on  the  sur> 
face  water  vapor  would  correot  the  integrated  refraction  computed  for 
dry  air  in  most  oases  but  would  almost  certainly  tend  to  give  divergent 
results  for  the  refraction  at  large  zenith  distances* 

The  contributions  to  the  total  refraction  by  individual  layers  for 
zenith  distances  of  60°  and  in  both  the  wet  and  the  dry  atmospheric 
model  are  given  in  Table  10.  For  z  -  60°,  the  layers  are  1  km  thick 
to  5  km,  and  5  kra  thick  from  IJ  to  30  km.  For  z  ■  85°>  the  layers  are 
1  km  thick  up  to  10  km,  and  5  km  thick  from  10  to  30  km.  The  contribu¬ 
tion  by  the  5  km  layer  between  50  and  55  km  is  also  shown  for  z  -  85°. 

The  differences  between  the  refractions  in  the  wet  and  the  dry  atmosphei*ic 
models,  acciunulated  down  to  any  level,  h,  is  shown  in  Figure  6.  It  is 
to  be  noted  that  from  their  effective  tops  (90  km)  down  to  30  km,  both 
atmospheres  give  essentially  the  same  amount  of  accumulated  refraction, 
and  the  maximum  differences  at  the  surface  are  only  0*l6U  in  the  z  ■  60° 
case  and  l!!o73  in  the  z  ■  85°  case.  The  total  refraction  confuted  by 
numerical  integration  is  given  in  Table  11  for  the  wet  and  the  dry  cases 
for  zenith  distances  of  60°  and  85°.  For  the  sake  of  ccnparison  the  values 
of  the  refraction  for  these  same  cases  confuted  from  the  Pulkova  tables 
and  using  Willis'  method  are  also  presented  in  this  table.  In  order  to 
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ooB|>uta  th«  values  for  the  wet  atnosphore  by  these  last  two  aethodSf 
the  surface  water  vapor  content  was  taken  to  be  that  assuasd  in  the 
•hiiolst  AROC"  ease  and  the  surface  tea^jerature  taken  to  be  that  required 
by  the  assuiqption  therein  of  an  unchanged  pressure  profile. 


Table  10 


REFRACTION  IN  INDIVIDUAL  LAIERS 
OF  THE  ARDC  (VIET)  AND  ARDC  (isr)  ATMOSPHERIC  H3I1BLS 


z  -  60°  »  -  85° 


Layer  (km) 

Wet 

Dry 

Layer  (km) 

i _ 

Wet 

Dry 

30  -  25 

l!l8l27 

i!!8127 

^5 

50 

0.^32 

o!!iU32 

25  -  20 

3.8527 

3.8526 

30 

- 

25 

8.I1608 

8.1*601* 

20  -  15 

8.U7U8 

8.U7U8 

25 

- 

20 

18.8939 

18.8933 

15  -  10 

17.5715 

17.57U 

20 

- 

15 

U3.9055 

1*3.9039 

10-5 

26.0232 

26.0161 

15 

- 

10 

96.72U7 

96.7188 

5  -  U 

6.6989 

6.692U 

10 

- 

9 

2U.6U51 

2U.6U32 

U  -  3 

7.2726 

7.2605 

9 

- 

8 

27.1*151 

27.1*098 

3  -  2 

7.881tU 

7:8623 

8 

- 

7 

30.1*31*0 

30.1*261* 

2-1 

8.53U7 

8.U9U3 

7 

- 

6 

33.7160 

33.7018 

1-0 

9.2390 

9.16U0 

6 

- 

5 

37.2736 

37.2525 

5 

h 

a.l597 

UI.II82 

U 

- 

3 

1*5.371*9 

1*5.2972 

3 

- 

2 

1*9.9676 

1*9.8252 

2 

- 

1 

5U.9589 

5U.6973 

1 

- 

0 

60.1*719 

59.9802 

I 
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Table  U 


SUMMARY  OF  REFRACTION  CALCULATIONS: 


Surface  Conditions 

Wet 

Qry 

t 

1U.19®C 

15°C 

p 

760  mm  Hg 

760  mm  Hg 

w 

.OOI16U 

0 

Method 

z  ■  60° 

ARDC  Integration  (A) 

98!l78U 

98!620 

*  0I16U 

Willis  Tables  (W) 

98.729 

98.563 

4-  0.166 

Pulkovo  Tables  (P) 

98.675 

98.512 

+  0.163 

A  -  W 

+  o!!o55 

♦  o!!o57 

A  -  P 

♦  0.109 

♦  0.108 

W  -  P 

♦  0.05U 

♦  0.051 

Method 

z 

00 

i 

ARDC  Integration  (A) 

579!!887 

578.81U 

♦  1I073 

mils  Tables  (W) 

580.176 

579.026 

♦  1.150 

Pulkovo  Tables  (P) 

579.657 

578.512 

♦  1.1U5 

A  -  W 

-  o!l289 

-  01212 

A  -  P 

♦  0.230 

0.302 

W  -  P 

+  0.519 

♦  O.51U 

The  agreement  between  the  results  of  the  thriie  methods  ie  ronark- 
ably  goodj  even  in  the  very  demanding  case  of  z  ■  85°  the  maxlniim  spread 

N 

between  all  three  methods  is  0.077  in  a  total  refraction  of  nearly 

580" .  The  85°  integration  Indicates  an  increase  of  the  wet  over  the 

a  „ 

dry  refraction  of  about  0.08  less  than  the  increase  indicated  by  the 

I 

other  two  methods,  but  the  actual  refraction  derived  from  the  integraticm 


falls  between  the 
cases. 


It  is,  of  course^'  important  to  note  that  all  values  calculated 
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are  baaed  upon  atmospheric  models}  one^  the  ARDC  atmosphere  being  an 
explicit  models  while  the  other  two  are  implicit  in  the  Pulkovo  and 
Willis  tables  for  conqputing  astronomical  refraction*  The  fact  that 
both  the-^ulkovo  and  Willis  tables  are  based  upon  precise  astronomical 
observations  and  that^  furthermore^  these  tables  can  be  used  in  the 
inverse  sense  to  coiipute  refraction  corrections  to  precise  astronomical 
measurements  implies  that  the  iitplicit  atmospheric  models  associated 
with  these  tables  approximate  exceptionally  well  the  actual  atmosphere^ 
at  least  on  the  average*  The  almost  negligible  divergence  of  the  re¬ 
fraction  integration  based  on  the  ARDC  atmosphere  from  the  Willis  or 
Pulkovo  values  implies^  at  the  leasts  that  the  ARDC  atmosphere  is  it¬ 
self  as  good  an  approximation  to  the  true  atmosphere  as  the  other  two* 
This  favorable  con^arison  also  implies  that  it  may  well  be  possible 
to  obtain  even  more  accurate  values  of  atmospheric  refraction  by  direct 
calculation^  that  iS|  integration^  of  an  atmospheric  model  representing 
a  suitable  modification  of  the  ARDC  models  the  modification  being  based 
upon  a  limited  sounding  of  the  atmosphere  at  the  observing  site  or 
(even  better)  upon  an  analysis  of  synoptic  weather  data  for  the  general 
locale  of  the  site* 
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VI.  RECOMMENDATIONS 

Th*r«  are  now  available  several  possible  ways  of  calculating  with 
considerable  precision  the  angular  refraction  in  a  vertical  plane 
caused  by  the  earth's  atmosphere.  These  range  from  the  astronoaioal- 
observation-based  Pulkova  Tables  used  in  one  or  another  of  their 
successively  ii^nrcved  forme  by  positional  astronomers  for  over  a 
century^  to  integration  through  the  idealised  ARDC  atmosphere  which 
has  been  derived  from  balloon,  rocket  and  satellite  soundings.  Viewed 
in  this  light,  the  really  remarkable  agreement  exhibited  by  the 
selected  eases  we  have  treeted  in  this  report  and  our  earlier  one,  is 
on  the  one  hand  assurance  that  any  one  of  the  several  diverse  approaches 
provides  accuracy  sufficient  for  the  vast  majority  of  present  needs, 
and  on  the  other  hand  a  strong  indication  that  we  are  indeed  also  in 
a  favorable  position  to  attenpt  some  further  gains  in  accuracy  through 
closer  adaptation  of  the  chosen  model  atmosphere  to  the  actual  condi¬ 
tions  obtaining  along  the  ray  path. 

1.  We  reecgmeend  that  a  careful  and  extensive  interoomparlson  be 
made  of  the  refraction  values  yielded  by  the  most  reomt  ARDC,  Tropical, 
and  Arctic  Winter  atmosphere  models  and  the  Pulkova,  Willis  and 
Astronomische-Cleodatische  Jahxbucb  (1952)  tables.  This  interooi^Mirlson 
will  seinre  two  objectives.  First,  it  will  provide  a  means  for  assessing 
the  magnitude  of  the  error  which  can  occur  in  the  deduced  position  of 
a  low  satellite  as  a  result  of  an  Incoiqpatiblllty  of  the  measured 
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avaettom  rcfraetiv*  iiid«x  vith  th«  true  total  atnoapharie  rafraetion 
at  tha  inatant,  an  errot  which  wa  hava  reason  to  baliara  la,  in 
practioa,  eonaidarably  lass  in  Magnitude  than  some  inrastigators  hare 
auggaatad.  Saoond,  it  will  provide  a  basis  for  the  study  of  rafraotlons 
in  the  "Iqrbrid*  aiodala  described  in  Raeomnendation  5* 

2,  Wa  raeonaMnd  that  tha  basic  Pulkova  tables  be  carefully 

sMoothed  to  one  more  daolMsl  plaea  and  tha  auxiliary  fomulas  and  * 

corresponding  tables  be  extended  by  one  decimal  through  utilisation  of 

tha  slightly  ii^troved  physical  constants  (a.g.  the  refractive  index 
of  water  vapor)  which  hava  becone  available  since  tha  nost  recant  edition 
of  tha  Pulkova  Tables.  This  will  reduce  the  number  of  "sensitive" 
cases  (for  low  satellites  or  very  large  zenith  angles)  where  round-off 
error  can  render  the  cospited  refraction  slightly  inconsistent  with 
the  surface  refractive  index  supposed  to  go  with  it  in  the  particular 
tabularly-defined  model.  A  similar  smoothing  (to  resMve  round-off 
error  but  not  to  affect  the  known  physical  discontinuities  incorporateid 
in  the  tabulated  models)  should  be  carried  out  on  the  AIOIC  and  related 
models  for  the  same  purpose,  to  the  extant  deemed  necessary  to  produce, 
for  the  sensitive  oases,  the  maxinun  accuracy  of  model  representation, 
attainable. 

3.  Ms  recoamwnd  that  an  integrated  study  be  made  of  the  systematic 
corrections  to  the  standard  astronomical  refraction  tables.  Such  correc¬ 
tions  have  already  been  derived  at  a  nudber  of  astronomioal  observatories. 

These  corrections  should  be  correlated  with  the  results  of  the  Interooqparisons 
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obtained  in  ReocoMndation  1.  These  correlations  should  then  aake  it 
possible  either  to  eetablieh  the  oxirrently  used  tables  (with  the  saall 
nodifioations  cited)  as  providing  adequate  accuracy  singly  or  in 
Judicious  coii>ination  for  all  presently  envisioned  neede,  or  to  derive 
(using  then  as  a  guide)  aethods  for  adjusting  chosen  physical  aodel 
atmospheres  singly  or  in  coidbination,  to  local  conditions,  along  the 
lines  su^geated  by  our  work  on  the  introdwtion  of  water  vapor  into 
the  ARDC. 

It*  We  recoanend  that,  along  with  this  study  involving  the  systeaetio 
"observatory  corrections,"  an  investigation  be  undertaken  of  selected 
portions  of  the  vast  aaount  of  data  on  the  more  or  less  random  varia¬ 
tions  in  refraction  available  for  the  groiq>  of  observatories  chosen 
under  Recoenendation  3*  Such  an  investigation  should  result  in  a 
meaningful  assessment  of  the  extent  of  the  scatter  of  refraction  effects 
in  actual  precise  observations  and  the  character  of  the  scatter  in 
relation,  for  example,  to  the  instantaneous  surface  index  of  refraction. 
This  in  turn  should  indicate  the  range  of  inoooqMtibility  of  that 
refractive  index  and  the  instantaneous  total  atmospheric  refraction 
actually  existent. 

$,  We  reconsmd  that  a  thorough  investigation  be  siada  of  the 
iaq>rovement  obtainable  in  "instantaneous"  totaH  ataospherio  refraction 
values  at  a  station  by  appropriately  coidaining  segments  of  the  ARDC, 
Tropical,  and  Arctic  Winter  physical  models,  suitably  adjusted  for  local 
variation  in  the  governing  paraateters.  These  atmoapheric  segments  should 
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b«  eoitbincd  by  t«ohnlqaM  alailar  to  thoa«  deaerlbed  in  tha  body  of 
this  report.  The  basis  for  the  selection  of  the  etaospherle  segaents 
would  be  the  air-issss  data  detendned  from  synoptic  weather  cherts  for 
the  surface  and  upper  air  relevant  to  that  particular  portion  of  the 
earth's  atmosphere  traversed  by  the  ray  being  studied*  Vfe  feel  that 
this  is  the  most  promising  mode  of  attack  for  obtaining  highly  precise, 
instiintaneous  values  of  the  refraction  which  are  truly  cosQiatible  with 
the  surface  parameters. 

6.  We  recommend  initiation  of  a  study  of  horisontal  refraction 

in  the  earth's  atmosphere  based  upon  both  the  observatory  data  mentioned 
above  and  the  geometry  of  the  asimuthal  discontinuities  at  air^sass 
fronts,  idiich  violate,  spatially  and  teaqporally,  the  assuaqpiion  of  a 
spherically  stratified  atsrasphere  basic  to  the  standard  foratulations 
for  vertical  refraction.  Such  a  study  could  be  pursued  most  efficiently 
as  an  adjunct  to  the  very  Important  one  outlined  in  Recommendation  5, 
since  it  would  Involve  much  of  the  same  information  and  some  of  the 
same  lines  of  thought. 

7.  We  recommend  that  the  effect  of  the  variation  of  refractive 
index  along  the  ray  path,  in  terms  of  both  geometrical  path-lei^h  and 
transit  tiam,  for  the  case  of  trllateratlon  on  a  satellite,  be  carefully 
studied  by  integration  methods  essentially  paralleling  those  we  have 
used.  The  three  atmospheres  (ARDC,  Tropical,  Arctic  Winter)  should  be 
taken  as  the  basic  models.  Two  problems  should  be  given  speoial  attentiont 
(1)  adaptation  of  the  atmospherio  model  to  take  account  of  the  air  masses 
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prMant  and  (2)  tha  aodlfioatlon  dua  to  introduction  of  watar  Taper. 

Tha  lattar  problea  la  particularly  iaportant  at  radio  fraquanelaa.  ro 
tha  flrat  ordar,  at  laaat,  such  a  study  is  stap-by-stap  an  analog  of 
tha  anrular  dariatlon  casa;  tha  intagratlon  would  prasuasbly  hara  to 
ba  axtandad  to  graatar  heights  because  of  tha  cunulatlTa  affect  on  tha 
transit  tlae,  of  tha  axtanslTa  regions  where  tha  index  is  not  quite 
unity  but  so  nearly  constant  that  no  aeasurabla  angular  affect  arises. 

8.  Wa  recosnend  that  careful  consideration  ba  glwan  to  the 
question  of  exactly  how  the  index  of  refraction  of  huodd  air  is  related 
to  tha  ware  length  of  tha  radiation  trawarsing  it.  In  tha  Tlslbla  and 
near  infra-red  the  index  decreases  slowly  with  Increasing  water  Taper 
content,  but  at  radio  waTe  lengths  tha  index  increases  rapidly  with 
Increasing  water  vapor.  It  is  therefore  possible  that  at  one  (or  more) 
intermediate  ware  length  the  effect  is  aero,  that  is,  that  the  refractiTe 
index  is  Independent  of  the  amount  of  water  vapor  present.  If  this 
should  be  the  case,  then  a  laser  (or  a  microwave  iMser)  operating  at 
that  wave  length  would  be  of  enoraious  potential  value. 


John  E.  Merrill,  Principal  Scientist 
kVitotamj 


Approved  byt 


id 


Francis  L.  Jacksm 
Director  of  Laboratories 
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